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Abstract 

We report classical and tight-binding molecular dynamics simulations of the Cgo fullerene and 
cubane molecular crystal in order to investigate intermolecular dynamics and polymerization pro- 
cesses. Our results show that, for 200 K and 400 K, cubane molecules remain basically fixed, 
presenting only thermal vibrations, while Cqo fullerenes show rotational motions. Fullerenes per- 
form "free" rotational motions at short times 1 ps), small amplitude hindered rotational motions 
(librations) at intermediate times, and rotational diffusive dynamics at long times 10 ps). The 
mechanisms underlying these dynamics are presented. Random copolymerization among cubanes 
and fullerenes were observed when temperature is increased, leading to the formation of a dis- 
ordered structure. Changes in the radial distribution function and electronic density of states 
indicate the coexistence of amorphous and crystalline phases. The different conformational phases 
that cubanes and fullerenes undergo during the copolymerization process are discussed. 

PACS numbers: 68.35.bp 71.20.Tx 82.20.Wt 87.10.Tf 
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I. INTRODUCTION 



The discovery of Ceo fullerenes^- opened a new field in theoretical and experimental re- 
search of carbon based materials^. At room temperature these molecules crystallize in a 
face-centered-cubic (fee) solid structure^ exhibiting semiconducting behavior (bandgap of 
about 1.5 eV-). When the crystals are either hydrostatically compressed^^^ or exposed 
to visible or ultra-violet light sources', a polymerization process of the Cqq molecules is 
observed. Due to the large size of the Cqo molecules, the interstitial cavities in the crys- 
talline Cqo can accommodate various guest species. Attempts were made to intercalate Ceo 
crystals with electron donors (and other atoms) in order to vary their electronic properties. 
Results of doping with alkaline metals lead to the appearance of exotic superconductivity 
with critical temperatures around 30 K or larger— Other atoms and molecules, such 
as rare gases" and molecular oxygen' were also used as dopants. 

Recently, heteromolecular crystals of Cqq fullerene and cubane (CgH 8 )' have been 
prepared from aromatic solutions by evaporating solvent and adding isopropyl alcohol 
as precipitant'. These crystals (Fig. 1) present an interesting phase diagram showing 
orientational-ordering phase transitions. At atmospheric pressures and temperatures be- 
low 140 K, the crystal exhibits an orthorhombic symmetry. From 140 K up to 470 K, the 
Cqq molecules are located on the lattice sites of an fee crystal with the cubane molecules on 
the octahedral voids. In this temperature range, the Cgo molecules are free to rotate whereas 
cubanes behave like static bearings in a so-called rotor-stator phased Above 470 K, solid- 
state reactions take place among cubanes and fullerenes leading to different crystalline and, 
eventually, amorphous phases due to copolymerization'. Similar phase diagrams are also 
observed for other fullerenes, such as C70 and Cs4— . 

Despite experimental studies about the polymerization process that takes place during 
high temperature treatments" and under pressure^, providing information of the possible 
conformations of the resulting polymer, the molecular aspects involved in these processes 
are still largely unknown. For instance, experimental results suggest that the resulting 
compound after thermal treatment (470 K) is a copolymer of Cqo with a decomposition 
product of cubane (dihydropentalene)'. However, other possibilities for the copolymer can 
not be ruled out even when other decomposition products of cubane are considered. In 
this case, similar copolymers could also be formed leading to different properties. While 
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detailed information on such processes is not easily accessible from experiments, atomistic 
simulations can be a very important tool to provide additional insights and molecular level 
understanding of the experimental data. In order to investigate the dynamical processes 
associated with Cgo-cubane crystals we present here atomistic simulations of the Cgo-cubane 
heteromolecular crystal. 

Firstly, we have analyzed the Cgo-cubane crystal in its rotor-stator phase. Our aim was 
to provide a detailed characterization of the Cqq and cubane intermolecular dynamics in 
the crystal. Secondly, we have studied the polymerization process in order to identify the 
possible conformations of the resulting copolymer and the changes that may eventually occur 
on the electronic properties of the crystal. 



II. METHODOLOGY 



The atomistic simulations were performed at two distinct levels: classical and tight- 
binding approximations. We have used classical molecular dynamics (MD) to investigate 
molecular motions in the rotor-stator phase. The use of such type of approximation is well 
justified because in that phase the main interaction between Cgo and cubane molecules is due 
to weak van der Waals forces. Furthermore, no covalent bondings among the molecules have 
been observed in the temperature range covered by the rotor-stator phase. The used force- 
field for the atom and bond interactions was the well-known Charmm parameterization^ 
in the Namd parallel molecular dynamics code^. In this type of simulation we investigated 
a supercell composed of 3x3x3 fee unit cells of the molecular crystal (108 fullerene and 108 
cubane molecules). Minimization of this system resulted in a lattice parameter of 14.50 A, 
only 1.6 % smaller than the experimental value (14.74 A—). The system was equilibrated 
during 50 ps in the canonical ensemble (NVT) followed by 200 ps runs in the microcanon- 
ical ensemble (NVE), where the molecular dynamics were analyzed. This equilibration run 
was sufficient for the system to reach thermodynamical equilibrium with the velocities fol- 
lowing the Maxwell-Boltzmann distribution. We considered two temperatures (200 K and 
400 K) within the range of the rotor-stator phase. Equilibration was performed carrying 
out Langevin dynamics with a damping coefficient 1 ps~ x with the Newtonian equations 
integrated with the Briinger-Brooks-Karplus method^. In the NVE simulations the Newto- 
nian equations of motion were numerically integrated by using the Verlet algorithm^ with 
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a timestep of 0.5 fs. 

The classical approximation using Charmm is no longer valid to study the polymerization 
process due to the existence of the possible bonding breaking/formation not possible in the 
Charmm approach. We have initially performed calculations with the well-known Tersoff- 
Brenner empirical potential designed for carbon systems^, which allows bonding breaking 
and formation during the evolution of the system, but it did not correctly describe the 
cubane geometry. So in order to investigate the polymerized phase we have described the 
interaction among atoms with a tight-binding model. This approach is intermediate between 
the classical empirical models and the fully ab initio density functional theory. Due to the 
size of the system investigated here for studying the polymerization process (304 atoms 
in the unit cell; C272 H 32 ) the tight-binding approximation presents the best compromise 
between accuracy and computational cost. 

We have then carried out tight-binding molecular dynamics (TBMD) simulations to in- 
vestigate the polymerization process at different temperatures. We used the tight-binding 
model developed by Porezag et al— implemented in the Trocadero program^. This 
model includes explicitly the non-orthogonality of the s-p basis, in which the hopping ma- 
trix elements are obtained directly from density functional theory calculations using the 
same basis set but disregarding three-center contributions to the Hamiltonian. This method 
has been successfully applied to the prediction of allotropic forms of carbon^ and polymer- 
ization of Ceo inside carbon nanotubes^ and it has been proven to combine accuracy and 
reduced computational effort, specially for large systems. 

Due to the high computational cost of TBMD, a single fee unit cell has been simulated 
during 40 ps (timestep of 0.5 fs) at 2100 K, 2300 K, 2600 K, and 3100 K. These temperatures 
are higher than those of the polymerization experiments, i.e, 470 K to 680 K— >2i Such high 
temperatures are needed in the TBMD simulations in order to accelerate the polymerization 
process (experimentally taking from seconds to hours for completion) to accessible timescales 
of the order of hundreds of picoseconds within current computational capabilities. Exper- 
imental results indicate that the polymerization reaction yields to only a slight expansion 
of the lattice^^i. Based on that, we considered the volume of the system fixed during the 
simulations using the lattice parameter obtained from minimization of a fee unit cell within 
our tight-binding model. The lattice parameter obtained was 15.08 A, only 2.3 % larger than 
the experimental one. We attribute this small difference to the not very accurate description 
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of van der Waals interactions within our tight-binding model. This is not a concern in the 
sense that for the high temperatures considered here, non-bonded interactions should play 
a minor role. The equations of motion were numerically integrated using the Generalized 
Leap- Frog algorithm^ and the temperature was controlled by means of the Nose-Poincare 
thermostat^. Only the T-point has been used for Brillouin-zone sampling. 

III. RESULTS AND DISCUSSIONS 
A. I. Rotor-stator phase 

The obtained results with classical MD simulations confirm the rotor-stator behavior of 
the molecular crystal, where the cubane molecules remain basically fixed (presenting only 
thermal vibrations) while Ceo fullerenes show rotational motions. As expected, we have not 
observed any translational motion of any molecule in the crystal. In Fig. 2 we present a 
mapping of the visited positions from a carbon atom of the Ceo and of the cubane during 
200 ps at 400 K. We can see that each atom of the cubane molecule only oscillates around 
its initial equilibrium position, whereas a targeted carbon atom of Ceo walks through all 
the allowed spherical surface for the same period of time. The moviel in the supplemental 
material provides a better visualization of these motions^. 

The observed movement for the fullerenes presents very interesting aspects. It is not a 
simple free rotation around a specific (and fixed) axis but a composed rotation instead, i.e., 
the rotation axis varies continuously with time. In this sense the rotor only executes a small 
fraction of a full cycle until the next change in the orientation of the rotational axis. In this 
case, due to environment presented by the neighboring cubane molecules, each fullerene atom 
develops (not independently) a random walk in angle over the allowed spherical surface. 

In order to better characterize the molecular dynamics of Ceo and cubane in the rotor- 
stator phase we computed the second-rank single-particle time correlation function C(r) 
defined as 



Here iV is the number of Ceo molecules, ^(x) is the Legendre polynomial of order 2, and 
u n (t) is the normalized orientational vector connecting the nth center of mass with an 
arbitrary atom in the nth Ceo at time t. The average is performed over different time origins 




(1) 



n=l 
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to- The C(t) time correlation function is useful to characterize reorientational motions 
in molecular systems and its Fourier transform is related to Raman- and light-scattering 
spectra^ 4 -. 

Fig. 3 depicts C(r) for the two temperatures considered here. For comparison, Fig. 
3 also shows C(r) computed analytically for a spherical free rotor— at 200 K with the 
Ceo moment of inertia. The reorientational time-correlations for fullerenes and cubanes are 
markedly different. The C(r) function for the cubanes exhibits rapid damped oscillations at 
short times (see insert) but no long time relaxation within the simulation time, indicating 
that the cubanes perform small amplitude fast librational motions (period ~ 0.2 ps) but 
no tumbling, consistent with the mapping pictured in Fig. 2. The fullerenes, in contrast, 
exhibit nearly full reorientational relaxation within this time span at both 200 and 400 
K. The fullerene reorientational relaxation in the Cgo-cubane crystal at these temperatures 
differs significantly from the free rotor behavior. The behaviors of fullerene and the free 
rotor coincide only at short times, up to about 0.6 ps (see insert in Fig. 3). 

Overall, the C(r) function reveals three distinct relaxation regimes for the fullerene 
molecules in the crystal, which resemble in many ways the reorientational relaxation of 
molecular liquids. A very short-time inertial dynamics, that lasts up to r ^ 1 ps, is as- 
sociated with the fullerenes "ballistic" or "free" rotational motions^ 4 - 1 ^. The subsequent 
oscillatory features at intermediate times (1 < r < 10 ps) stem from small amplitude hin- 
dered rotational motions (librations) of the fullerenes under the influence of restoring forces 
or torques from the environment, largely from neighboring cubane molecules. The post- 
librational, long-time behavior (r > 10 ps) of C(r) is well described by an exponential 
decay characterizing a regime of rotational diffusive dynamics. Single exponential fits to 
C(t) at long times yield characteristic reorientation relaxation times of about 40 ps at 200 
K and 13 ps at 400 K for the fullerene molecules in these crystals. The overall relaxation 
time, r J , obtained from the time-integral of C(r), r J = J °° C(r)dr, which is related to 
the nuclear magnetic resonance (NMR) relaxation rate under extreme narrowing conditions, 
turns out 34 and 9 ps at 200 and 400 K, respectively. Our simulations suggest that the 
reorientational relaxation of the fullerenes is nearly one order of magnitude slower than that 
of a low molecular-weight nonpolar aromatic molecular crystal such as benzene^ (moment of 
inertia, I ~ 14xl0~ 46 kg.m 2 ) in this temperature range (r J = 0.5 — 6.0 ps), but somewhat 
more similar to that of carboranes (I ~ 5xl0 -45 kg.m 2 ), despite the differences between the 
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molecular rotational motions in these crystals 3 ^ 1 ^. 

It is also interesting to compare the behavior of Ceo in the Ceo-cubane crystal with the 
dynamics of the crystalline Ceo- Johnson et al. have investigated the rotational dynamics 
of Cgo in the solid state using NMR over the temperature range of 240 up to 331 K— . 
They observed that the reorientational correlation time follows an Arrhenius temperature 
dependence with activation energy of approximately 1.4 kcal/mol for the rotor phase of 
crystalline Ceo- The value of the activation energy is of the same order of the energy barriers 
found in the rotational potential maps (Fig. 5). If we consider the rotor phase in the fee 
structure (temperatures above 260 K— ) and extrapolate the NMR temperature dependence 
line of the reorientational correlation time in^ up to a temperature of 400 K, we obtain 4.6 
ps^i. The corresponding correlation time for the free gas Cgo is about 2.5 ps^. Therefore, 
our simulations indicate that the reorientational relaxation time for Ceo-cubane crystal is 
about two times as long as the expected time for the Ceo hi the solid state, and almost four 
times as long as the estimated time for unhindered gas-phase rotation at around 400 K. 

Fig. 4 depicts the frequency spectrum of C(r) for the two temperatures investigated here. 
Two peaks are evident in the spectra: one at 6 cm -1 and another at 14.5 cm" 1 . The latter 
is more evident in the case of 200 K while the former is more pronounced at 400 K. We 
associate these peaks mainly with the molecular reorientation within the librational regime. 

In order to obtain insights into the origin of such peaks we first calculated the rotational 
potential energy map for the fullerene due to van der Waals bare interactions with its six 
neighboring cubane molecules (Fig. 1). For these calculations we kept the cubane molecules 
fixed and varied the angular orientation of the fullerene, neglecting any libration-phonon 
coupling. The angular orientation of the fullerene inside the cavity (Fig. 1) is characterized 
by the Euler angles (0 < < 27T, < 9 < n, and < tp < 2ir), and determined by 
the rotation matrix R(0, ip). For each combination ((f), 6, ip) the potential energy U was 
calculated leading to a four-dimensional map (U (cp, 9, ip)). The three-dimensional potential 
energy map calculated for each value of <p was then used to determine local minima. For 
each minimum we calculated the corresponding rotational frequency u r through 

1 d 2 U 



,2 



(2) 



21 d£ 2 

where / is the Ceo moment of inertia (9.75 x 10~ 44 kg.m 2 ) and £ = tp, 9. These minima 
represent local rotational traps for the fullerenes in the crystal. We can see from Fig. 5 
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that different minimum valleys appear depending upon the fullerene orientation. Fig. 6(a) 
shows the uj r distribution obtained from the potential energy mapping of a Ceo fullerene 
surrounded by six cubane molecules (Fig. 1) for 24 different values of <fi. The distribution 
shows frequencies up to about 10 cm -1 with a maximum around 6 cm -1 . The frequency at 
the maximum of the distribution corresponds to the frequency of the main spectral peaks 
in Fig. 4, thus showing compelling evidence that the main peaks of the reorient at ional 
frequency spectra (Fig. 4) are related to (hindered) rotational motions of the fullerenes 
under the potential energy wells due to the neighboring cubane molecules. We have not 
found any evidence of uj r values around 14 cm -1 in our limited search (24 cj) values). However, 
this simplified analysis neglects any effects due to lattice phonons (and internal vibrations 
of fullerenes and cubanes). It is physically sensible to expect that lattice phonons act as 
sources of torques thus influencing the librational motions of the fullerenes. In order to 
examine to what extend the librational motions of the fullerenes and lattice vibrations may 
be coupled, we have computed the time correlation function C„(t) for the fullerene center 
of mass velocity v cm (t): 

1 N 

C v (r) = — (v cm (to + r) ■ v cm (t )])i 0! (3) 

n=l 

and determined its spectral density C v (uS) using Fourier transformation. Results depicted 
in Fig. 6(b) show that there are indeed lattice vibrational modes in the 12-17 cm -1 region, 
thus supporting our interpretation. 

B. II. Polymerized phases 

Having described the dynamical behavior of the rotor-stator phase, we also investigated 
the polymerized phases. Cubanes can thermally isomerize into different compounds. In the 
experiments reported by Kovats et alM they showed that cubane molecules inside the Cqo- 
cubane crystals isomerize at the same temperature (about 200°C) as in their free standing 
forms. Kovats et al. also suggested that, among the most frequent decomposition products 
of cubane^i^ (Fig. 7), dihydropentalene (DHP) and styrene (STY) appear as the most prob- 
able compounds to form alternating copolymers with Cgo- Table I presents the energy gain 
for the cubane isomerization into different products obtained from pyrolysis experiments^ 
and also from the present calculations. We can see from the Table I that STY appears 
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as the most stable product, followed by DHP. This is also predicted by the tight-binding 
model used here, as well as, by more sophisticated density functional theory calculations. 
Furthermore, both types of calculations predict the same stability ordering. The calcula- 
tions are consistent with the experimental values, except for the bicyclooctatriene (BCT) 
and cyclooctatetraene (COT) products which appear exchanged. 

Martin et al. experimentally investigated cubane in the temperature range 230— 260°C 
having observed its decomposition to COT and further fragmentation to benzene (BEN) 
+ acetylene (C2H2) at low pressures^&^i. Experimental results also indicate that fur- 
ther products of COT are DHP at 700-850 K and STY and BEN + C 2 H 2 at higher 
temperatures^. These results suggest that despite STY and DHP being the most stable 
products, the isomerization pathway has COT as an intermediate product in these pro- 
cesses. 

In order to analyze the cubane decomposition process in its gas phase, we have carried 
out tight-binding molecular dynamics simulations of a system composed by eight cubane 
molecules during 50 ps at different temperatures. These simulations also allowed us to 
estimate, at least approximately, the simulation temperature which cubane molecules begin 
to isomerize. In the simulation time scale used here, we have verified that at 1000 K no 
isomerization occurs. When the temperature is increased to 1300 K, COT molecules are 
present, and at 1900 K COT eventually decomposes into BCT or DHP molecules. This 
process is consistent with experimental findings^ 1 ^^ 1 ^ and may be used to map simulation 
and actual experimental temperature scales. Thus, a temperature of ~ 1900 K in the 
simulations can be mapped into the actual polymerization temperature of ~200° C. We 
analyzed the Cgo-cubane crystal at temperatures larger than 1900 K which ensures, therefore, 
the decomposition of cubane molecules. 

In fact, STCO, COT, BCT, and BEN + C 2 H 2 were already present at 2100 K after 10 
ps of TBMD simulations. Bonding between cubane products and Cgo begins to occur after 
20 ps, but they are not very stable. After 40 ps, the presence of COT is still observed as 
well as bonding among products and Ceo- Eventually bonding among Ceo molecules through 
products of the cubane decomposition were observed, as shown in Fig. 8 (a). 

After 10 ps at 2300 K all the cubane molecules have been decomposed and DHP com- 
pounds begin to appear. In addition, we have observed hydrogen atoms bonded to the Cgo 
molecules. After 40 ps, covalent bonding between BEN and Ceo as well as between DHP 
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and Ceo were observed as shown in Fig. 8 (b). When the temperature is further increased to 
2600 K, all the cubane molecules were decomposed into BEN + acetylene after 20 ps. The 
copolymerization of C 60 via cubane products is more easily observed after 30 ps as shown 
in Fig. 9 (a). 

At 3100 K, the extension of the copolymerization is much larger than that observed at 
2100 K but clearly shows signs of increased disorder. This can be seen in Fig. 9 (b) where 
three fullerenes are connected through fragments of cubane molecules. In this case, the 
crystal shows features reminiscent of an amorphous phase. We also observed that some 
fullerene cages were opened during the polymerization process at 3100 K. 

Based on these results, crystalline and amorphous aspects are present for 2100-2600 K 
while a solid with a more pronounced amorphous characteristic is observed for 3100 K. 
These findings agree with the experiments performed by Pekker et al.—. The evolution of 
the polymerization process with the temperature can be seen in Fig. 10 through the radial 
distribution function g(r) for 2100 K and 3100 K after 40 ps. For comparison purposes we 
also show g(r) for the crystal structure at K. The peaks at 1.39, 1.44, 2.46, and 2.85 A 
correspond to the first, and second nearest neighbors in Ceo, in agreement with neutron- 
diffraction measurements at 300 K 53 . The peak at 3.56 A corresponds to third nearest 
neighbor in Ceo, and at 1.12 and 2.35 A to the C-H and C-C bonds in the cubane molecules, 
respectively. 

The broadening of the main peaks is observed at 2100 K and, more significantly, at 
3100 K after 40 ps of molecular dynamics simulation. This indicates changes in the crystal 
structure due to bonding among Ceo and cubane products (2100 K) as well as bonding among 
C 6 o themselves (3100 K). The almost complete disappearance of the peak corresponding 
to the C-H bonds altogether with the broadening of peaks at 3100 K clearly show the full 
decomposition of the cubane molecules and the further polymerization of their decomposition 
products with the fullerenes. 

Finally, we present the observed changes in the electronic density of states (DOS) of 
the Ceo-cubane crystal during polymerization. Fig. 11 shows the DOS for three different 
temperatures. These results were confirmed with density functional theory calculations^. 
At K the crystal exhibits a semiconducting behavior with a bandgap of about 1.9 eV, 
where the electronic states near the Fermi level are mainly due to Ceo states. The DOS is 
significantly different at 2100 K and 3100 K. In both cases the DOS is broader than the 
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one obtained at K. Furthermore, an increased concentration of states at the Fermi level 
is clearly observed. These aspects correspond to signatures of the amorphous phase and 
disordered characteristics presented by the resulting structure after thermal heating. The 
band gap of 1.9 eV can be associated to a wavelength with 654 nm and consequently to an 
emission/absorption of an orangish light, consistent with experimental findings of Pekker et 
al. in Ceo-cubane crystals without heat treatment^. In addition, the broadening observed in 
the calculated DOS for 2100 K and 3100 K can indicate changes in the absorption properties 
of the crystals and might explain the experimental observations of darker crystal colors after 
heat treatment^. 



IV. SUMMARY AND CONCLUSIONS 



Classical and tight-binding molecular dynamics simulations were used to investigate the 
rotational dynamics and the polymerization processes in the Ceo-cubane crystal as function 
of temperature. Our results show that this crystal behaves like a rotor-stator system in 
agreement with experimental results^. The cubane molecules work as stators and fullerenes 
as rotors. The fullerene molecules perform free rotations only during small time periods. 
For larger time intervals small amplitude hindered rotational motions and random rotations 
were predicted. The characteristic reorient at ional time of Ceo m the Ceo-cubane crystal was 
predicted to be significantly larger than that exhibited by the Ceo crystal and by a fully 
free rotor. The fast components of the reorientational dynamics are characterized by two 
well- resolved bands in the frequency spectrum centered at 6 cm -1 and 14.5 cm -1 , which are 
associated with the curvature of the potential energy surface created by cubane molecules 
surrounding Ceo an d the coupling between Ceo hindered rotations and lattice vibrations, re- 
spectively. Random copolymerization is observed when the temperature is sufficiently high 
to allow cubane decomposition. The products of this decomposition initiate the polymeriza- 
tion process covalently connecting neighboring fullerenes. In agreement with experimental 
results^, dihydropentalene appears as an important cubane product in the polymerization 
process. However, we have also observed that cyclooctatetraene and benzene+acetylene 
(after styrene decomposition) also contribute to the copolymerization. The polymerization 
causes disorder in the crystal changing the local bonding environment and the electronic 
structure. The density of electronic states is significantly broadened and accompanied by 
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an increase in the concentration of states in the vicinity of the Fermi level upon copolymer- 
ization at high temperatures. We hope the present results help to interpret some unclear 
experimental data as well as on the design of new experiments to these structures. 
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TABLES 

TABLE I: Energy gain (in eV) of the cubane products obtained from experiments (EXP)^ and from 
our tight-binding (TB)2I and density functional theory (DFT)^ calculations. The experimental 
data were obtained from the differences in the heat of formation of the compounds. 



Cubane product 


EXP 


TB 


DFT 


Styrene 


4.92 


4.50 


2.79 


Dihydropentalene 


3.80 


3.24 


2.06 


Bicyclooctatriene 


3.02 


2.42 


1.33 


Cyclooctatetraene 


3.36 


0.65 


1.12 


Syn-tricyclooctadiene 


1.34 


0.23 


0.29 
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FIGURE CAPTIONS 

Fig.l: Hetero molecular crystal structure of Ceo-cubane^. Each Ceo in the crystal is 
surrounded by six cubane molecules. In such molecular arrangement the cubane molecules 
work as static bearings and the Ceo behaves as a nearly free-rotating spheroid. 

Fig.2: (Color online) Mapping of the visited positions (smaller (blue) spheres) from a 
carbon atom of Ceo an d cubane during a 200 ps simulation at 400 K. The initial structures 
of the Cgo and cubane are superimposed to the mapping in order to help the visualization. 

Fig.3: (Color online) Single-particle correlation function C(r) vs r for the Ceo fullerene 
and cubane. For comparison purposes the C(r) for a free fullerene spherical rotor at 200 
K is also presented. The insert shows details of C(r) short-dynamics in a semi- log scale. 
The correlation function for cubane in the inset graph was expanded five times to facilitate 
visualization. 

Fig.4: u 2 times the frequency spectrum of C(r) obtained from the cosine Fourier 
transform (C(u)). As depicted, the frequency spectra would correspond roughly to the 
intermolecular components of the Raman or depolarized light scattering spectra of fullerenes 
in the crystal. 

Fig.5: (Color online) Rotational potential experienced by a Ceo fullerene when sur- 
rounded by six cubane molecules (Fig. 1) as a function of the 9 and ip angles for (a) = 0, 
(b) 30°, (c) 90°, and (d) 150°. 

Fig.6: (a) Distribution of the rotational frequency u r obtained from the curvature of 
the minima on the three-dimensional potential energy maps generated by rotating a Ceo 
fullerene surrounded by six cubane molecules around a fixed center of mass position (cf. 
Fig. 1). (b) Low frequency portions of the spectral density C v {u) of C v {t) obtained from 
Fourier transforms. 

Fig. 7: (a) Cubane and its most frequent decomposition products: (b) syn- 
tricyclooctadiene (STCO), (c) cyclooctatetraene (COT), (d) dihydropentalene (DHP), (e) 



1(3 



bicyclooctatriene (BCT), and (f) styrene (STY). 



Fig.8: Snapshots of the tight-binding molecular dynamics simulations showing fragments 
of the Ceo-cubane crystal at (a) 2100 K and; (b) 2300 K, after 40 ps. 

Fig.9: Snapshots of the tight-binding molecular dynamics simulations showing fragments 
of the Ceo-cubane crystal at (a) 2600 K and; (b) 3100 K, after 30 ps. 

Fig. 10: (Color online) Radial distribution function g(r) for the Ceo-cubane crystal at 
different temperatures. 

Fig.ll:Electronic density of states (DOS) for three different configurations of the Ceo- 
cubane crystal. The vertical dashed line indicates the position of the Fermi level. Non-zero 
density of states at the Fermi level for K is due to peak broadening. 
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